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We study the distribution of partition function zeroes for the Xy-model in two dimensions. In particular we 
find the scaling behaviour of the end of the distribution of zeroes in the complex external magnetic field plane in 
the thermodynamic limit (the Yang-Lee edge) and the form for the density of these zeroes. Assuming that finite- 
size scaling holds, we show that there have to exist logarithmic corrections to the leading scaling behaviour of 
thermodynamic quantities in this model. These logarithmic corrections are also manifest in the finite-size scaling 
formulae and we identify them numerically. The method presented here can be used to check the compatibility 
of scaling behaviour of odd and even thermodynamic functions in other models too. 



1. KOSTERLITZ THOULESS SCALING 

The partition function for the d-dimensional 
0{n) non-linear cr-model is 
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where L denotes the linear extent of the lattice 
A, /3 is the Boltzmann factor and h is the re- 
duced external field. The n-component spin at 
has unit modulus. In the case oi d — n = 2 this is 
the plane rotator or AF-model and has a phase 
transition at (/3, h) = {f3c, 0) which is caused by 
the binding and unbinding of vortices. 

An approximate renormalization group ap- 
proach indicates unusual scaling behaviour of 
the thermodynamic functions as criticality is ap- 
proached from the disordered phase m . In terms 
of the reduced temperature i = 1 — /3/f3c, the 
scaling behaviour of the correlation length, sus- 
ceptibility and the specific heat is given in H] as 
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where for i > 0, i^ = 1/2, rjc = 1/4 and a — 
—d = —2. The aim of this work is to argue that 
the latter two scaling formulae are incompatible 
as they stand. To this end, a method is presented 
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by which odd and even thermodynamic functions 
(hke (H) and ^) can be related and expressed 
in terms of partition function zeroes. Using cer- 
tain reasonable assumptions regarding finite-size 
scaling, it is shown that there have to exist mul- 
tiplicative logarithmic corrections to (y) and (||). 
This method can be applied to any model. 

2. LEE YANG ZEROES 

For the Ary2^inodel the zeroes in the complex 
external field strength plane are on the imaginary 
axis |2,g. Expressing the partition function in 
terms of its Lee- Yang zeroes (denoted Zj{f3)), 



ZL{l3,h) = pL{l3,h)l[{h - tz.m 
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where pL is a non-vanishing function of h related 
to the spectral density and contributes only to 
the regular part of the free energy. The singular 
part of the free energy corresponds to 

.fL{(3,h) = L-''Y.\n{h-zz,m . (6) 

j 
This can be written as 

/l(/3,/i)=/ lii{h~iz)gL{P,z)dz , (7) 

in which R is an appropriate cutoff and the den- 
sity of zeroes, 5l(/3, z) is given by 
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Here Gl is the cumulative density of zeroes. The 
distribution of zeroes is symmetric in the real h- 
axis and the (cumulative) density of zeroes is zero 
up to the Yang-Lee edge. This allows one to write 
the singular part of the free energy in the ther- 
modynamic limit as 
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where zi(/3) is the position of the Yang-Lee edge. 
This leads to an expression for Xoo in terms of 
the cumulative density of zeroes and the edge. 
Following [QJ5) this gives (independent of the 
model under consideration provided it obeys the 
Lee- Yang theorem) 
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where $ is unknown. The specific heat can also 
be written in terms of zi and Goo Wtn '^ 
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3. SCALING AND CORRECTIONS 



Assume the following modified Kosterlitz- 
Thouless (KT) scaling behaviour for the singular 
parts of the thermodynamic functions: 

xoo{t) - e^'^^f , (12) 

Coo(i) ^ Ct' ■ (13) 

Assume, furthermore, the edge scales as 
zi(t)~eitP • (14) 

Putting these in (nw and ^Wj gives M 

A = -(5-2 + 77,), P^^iq-r) + l + i^, (15) 

Thus the scaling behaviour of the Yang-Lee edge 
has multiplicative logarithmic corrections even if 
Xoo and Coo have not. 

4. FINITE SIZE SCALING 

Finite-size scaling (FSS) allows one to find 
the volume dependency of thermodynamic quan- 



tities at Pc from their thermodynamic limit scal- 
ing behaviour. Its general form (valid in all di- 
mensions including the upper critical one) is H 
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where Phit) is some thermodynamic function at 
reduced temperature t for a system of extent L 
and Fp is unknown. For the XY2-mode\ Cl(0) oc 
L §. Applying (|l|) to (|l|) and (0) gives 

Xl{0) ^ L^-^^{\nL)-i , (17) 

zi{L) - Li("-2+'k)(lnL)-i(^+i+'^) .(18) 



The finite-volume counterpart of (|10[), relating 
the susceptibility and the zeroes, is (from 
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where it has been assumed that the lowest lying 
zero has the dominant effect. This gives 



(5 = -d=-2, q = -2{1 + ly) = -3. 
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Thus the scaling behaviour of the singular part 
of the specific heat indeed exhibits multiplicative 
logarithmic corrections. 

5. NUMERICAL RESULTS 

The above arguments have yielded no informa- 
tion on the odd correction exponent r. Careful 
renormalization group (RG) considerations give 
r = —1/16 |lfl|. Our task now is to identify r 
numerically. To this end we use the FSS of the 
Lee- Yang zeroes. Accepting the KT predictions 
for v, rjc and a for i > 0, (|l8|) gives 
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A Wolff algorithm |llj was used to simulate 
the AF-model on square lattices of sizes L = 
32, 64, 128 and 256. The critical /3-value was 
found to be 1.11(1) from phenomonological RG 
methods [Q. The results for the lowest lying 
zeroes at (3 = 1.11 are zi(L) = 0.0023348(7), 
0.0006350(2), 0.00017279(5) and 0.000047062(13) 
for L = 32, 64, 128 and 256. Details on the nu- 
merics and the determination of f3c and the zeroes 
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Figure 1. Leading FSS of Lee- Yang Zeroes. 
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Figure 2. Corrections to FSS of Lee- Yang Zeroes. 



will be given in a forthcoming publication H . In 
the absence of any corrections, the slope of Fig- 
ure P should give the leading power-law FSS ex- 
ponent. In fact the slope is -1.8777(4), the devia- 
tion from the KT value of —15/8 = —1.875 being 
due to the presence of logarithmic corrections. To 
identify these, and the correction exponent r in 
(pT|), In (ziL^^/^) is plotted against Inlni in Fig- 
ure 0. A straight line is identified. Its slope is 
—0.012(1). Thus we have strong evidence for a 
non-zero value of r, albeit not in agreement with 
the RG predictions of -1/16 = -0.0625 from [|lO|. 



6. CONCLUSIONS 

Theoretical arguments checking the consis- 
tency of the scaling behaviour of odd and even 
thermodynamic functions at a KT phase tran- 
sition have been presented. The generally used 
scaling formulae have to be modified by multi- 
plicative corrections. These are identified ana- 
lytically for the specific heat and numerically for 
the susceptibility. This numerical identification 
comes via an analysis of Lee- Yang zeroes, the FSS 
of which is linked to that of the susceptibility. 

We would like to thank P. Lacock for assistance 
with multihistogram reweighting techniques. 
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